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Abstract 



O . We study the underpotential deposition of Cu on single-crystal Au(lll) elec- 

O 

trodes in sulfate-containing electrolytes by a combination of computational 
statistical-mechanics based lattice-gas modeling and experiments. The exper- 
imental methods are in situ cyclic voltammetry and coulometry and ex situ 
Auger electron spectroscopy and low-energy electron diffraction. The exper- 
imentally obtained voltammetric current and charge densities and adsorbate 
coverages are compared with the predictions of a two-component lattice-gas 
model for the coadsorption of Cu and sulfate. This model includes effective, 
lateral interactions out to fourth-nearest neighbors. Using group-theoretical 
ground-state calculations and Monte Carlo simulations, we estimate effective 
electrovalences and lateral adsorbate-adsorbate interactions so as to obtain 
overall agreement with experiments, including both our own and those of 
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other groups. In agreement with earher work, we find a mixed (\/3x\/3) 
phase consisting of 2/3 monolayer Cu and 1/3 monolayer sulfate at inter- 
mediate electrode potentials, delimited by phase transitions at both higher 
and lower potentials. Our approach provides estimates of the effective elec- 
trovalences and lateral interaction energies, which cannot yet be calculated 
by first-principles methods. 

1995 PACS Numbers: 68.45.-v 82.20.Wt 82.45.-Fz 81.15.Pq 
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I. INTRODUCTION 



Underpotential deposition (UPD) is the process by which a submonolayer or monolayer 
of one metal is electrochemically adsorbed onto another at electrode potentials that are 
positive with respect to the potential for bulk deposition. It occurs when the adsorbate 
adatoms are more strongly bound to the foreign substrate than to a substrate of their own 
kind This phenomenon has been known since the early years of this century 
but its study has intensified only during the last two decades [0,0. Current-potential data 
for UPD on single- crystal surfaces, obtained from cyclic voltammetry (CV) experiments, 
frequently display sharp current peaks 0J^. The shapes, positions, and number of these 
peaks depend on the substrate and the crystal plane on which the adsorption takes place 
0-0, as well as on the nature of the electrolyte ||lOl- It was already suggested in Refs. |]7|,|8|JTT[| 
that the CV current peaks separate ordered and disordered adsorbate structures, and it was 
later pointed out that the peaks themselves should correspond to phase transitions in the 
adsorbed layer |]l],|12|,|13[ . The UPD phenomenon is of fundamental interest as an aspect 
of the detailed, microscopic structure of electrode-electrolyte interfaces, which is presently 
the focus of vigorous research It provides a means of controlling the microscopic 

surface structure through the electrolyte composition and the applied potential, and likely 
applications are suggested by the observations that metallic submonolayers have the ability 
to significantly modify the catalytic activity of a substrate ||15|, and that UPD is the first 
step in the liquid-phase epitaxial growth of metallic heterostructures . 

The phase-transition aspects of UPD make this area of surface electrochemistry well 
suited for statistical- mechanical modeling 0-|2^. The first detailed, statistical- mechanical 
lattice-gas model for a UPD system was introduced by Huckaby and Blum (HB), who con- 
sidered the electrosorption of copper (Cu) onto the (111) plane of gold [Au(lll)] in the 
presence of sulfate ||T^ , p!3|pl| -p6| . This system also has been intensively studied by a number 
of experimental techniques (see Sec. |I|), and the detailed data that are now available render 
it an excellent candidate for theoretical modeling. 
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In this article we present a study of the UPD of Cu on well-characterized, single-crystal 
Au(lll) electrodes in sulfate-containing electrolytes, in which we employ both experimental 
and theoretical methods. The experimental techniques used are in situ CV and coulometry 
to detect the current peaks and measure the total charge transferred during the adsorp- 
tion/desorption processes, and ex situ Auger electron spectroscopy (AES) and low-energy 
electron diffraction (LEED) to determine the adsorbed Cu and sulfate coverages and the 
adlayer structure, respectively. In the theoretical part of the study we formulate a full, 
two-component lattice-gas model with effective, lateral interactions out to fourth-nearest 



neighbors. This is a generalization of the HB model [p^ , p!^ ,PH26| . Using group-theoretical 



ground-state calculations at zero temperature p7|-p9|l and Monte Carlo (MC) simulations at 



room temperature [p0|-|3^ , we estimate the model parameters (the effective, lateral interac- 
tions and electrovalences) so as to obtain overall agreement with the available experimental 
data, including both our own and those of other groups. This theoretical approach has 
previously been applied to electrochemical problems by Rikvold and co-workers p3|-p7|, and 
it was most recently used to model the coadsorption of hydrogen with urea on single- crystal 
platinum (100) electrodes [|T^,^^-^ and with sulfate on rhodium (111) JT^. It differs from 



the approach taken by HB in that we apply the full, two-component lattice-gas model to 
analyse global, as well as local features of the system's phase diagram and CV current 
profiles, and in that our room-temperature results are obtained by nonperturbative, numer- 
ical simulation. (For reviews and discussion of the HB approach, see Refs. p5|j26[] .) Some 
preliminary results of this study have been reported in Refs. pO| , |4^ . 

The remainder of this paper is organized as follows. In Sec. || we give a survey of 



experimental results that are currently available. In Sec. |T| we describe our experimental 
procedures and results. In Sec. ^ we introduce our lattice-gas model and give the relations 
between the lattice-gas parameters and the experimentally observable quantities. Here we 
also discuss our procedure for adjusting the model parameters to fit the experiments. In 
Sec. |V] we present together our numerical and experimental results, and in Sec. |VI| we give 



a finite-size scaling analysis of the phase transitions predicted by our model. In Sec. [VI 1| we 



briefly summarize our conclusions and discuss their significance. 



II. SURVEY OF THE EXPERIMENTAL SITUATION 

The UPD of Cu on Au(lll) in sulfate-containing electrolytes has been extensively 
studied, both by in situ techniques like CV [[7|,p|, p!0| , pA| , ^5| -pl[] , chronocoulometry ||i5|-|5T|, 



scanning tunneling microscopy (STM) |]9| j45|j5^j53[] , atomic- force microscopy (AFM) |^4| 
Fourier-transform infrared spectroscopy (FTIR) |^5|, X-ray-absorption spectroscopy (EX- 



AFS and XANES) [pq -|58|, surface X-ray scattering ||59|j60|l , quartz crystal microbalance 
(QCM) p7||59| , and by ex situ techniques, such as LEED, RHEED (refiection high-energy 
electron diffraction) and AES [|lO| , |rTl , |6ll -|63|] . Electrolyte compositions used in several CV 
experiments (including the present work) are listed in Table |, together with potential scan 
rates and observed current-peak separations. The observation that the different electrolyte 
compositions included in Table | do not qualitatively change the shape of the voltammogram 
indicates that within this range of compositions the adsorbates remain the same. In this 
paper we assume these to be Cu^"*" and S04~ in solution, which are adsorbed as near-neutral 
Cu and partially discharged or neutralized S0|^, respectively |^8|- ^|64 . 



A three-step experiment was performed by Shi and Lipkowski [|^,^ to clarify the roles 
of the different adsorbates. With 0.1 M perchloric acid (HCIO4) as the only electrolyte 
{i.e., with no copper added), the CV current in the double-layer range was weak, indicating 
that the ions present in the solution did not adsorb significantly on the surface. In the 
next step, K2SO4 was added, introducing HSO4 and S04~ ions into the solution. An 
adsorption peak then appeared, particularly during the positive-going scan. Radiochemical 
studies indicate that this peak is due to sulfate adsorption with a maximum coverage of 



approximately 0.2 monolayers (ML). This result agrees with chronocoulometric results 
and with earlier STM studies p5|j66|] , in which a (v^x^/7) sulfate adlayer was indicated 
on the positive side of the voltammetric peak in the absence of copper, without pH] or 



with a perchloric-acid supporting electrolyte. (In UHV studies with samples emersed 



5 



at potentials for which the (a/3x-\/7) structure is observed in situ, only a (-\/3xa/3) surface 
structure is observed, presumably produced by desorption of water that was coadsorbed 
with the sulfate [0].) Finally, addition of Cu^+ [as Cu(C104)2] caused the appearance of 
two sharp (10-30 mV wide) current peaks that totally dominate the voltammogram, as 
illustrated by our own data in Fig. |l|(a). There is no two-peak feature in the CV if no 
copper is added. On the other hand, if the copper concentration is kept fixed while the 
concentration of sulfate is reduced, either the peak heights are reduced, or the two peaks 
merge, or both effects are observed together |^ , |49[| . Evidently, the coadsorption of Cu and 
sulfate is mutually enhancing. (The same effects, albeit with a smaller peak separation, are 



also seen with Pt(lll) electrodes in the same electrolytes ^8 



In Fig. |I](b) is shown the charge density transferred during the adsorption/desorption 
process, obtained as the integral of the CV current density in Fig. |I|(a) between +120 and 
+420 mV vs. Ag/AgCl. The charge plateau in the potential region between the CV peaks 
(approximately 200-300 mV) lies about 200 fiC/cm^ below the charge density at 420 mV. 
If the reactions Cu — 2e^ Cu^+ and SO4 + 2e^ S0|^ are assumed across the right- 
hand peak (Peak 7^1), the desorption of 2/3 ML of Cu and 1/3 ML of sulfate would result 
in a total charge transfer in that potential region of 160 /iC/cm^ [^. This coverage- 
charge correspondence is supported by RHEED and LEED llT],^, by which a (a/3 x a/3) 
superstructure was observed in the potential region between the peaks. Further, consistent 
results have been obtained by STM Corroborating evidence was also obtained 

in an in situ electrochemical QCM study At the left-hand peak (Peak #2) and for 

still more negative potentials, the Cu coverage approaches 1 ML, and EXAFS results have 
suggested that sulfate adsorbs on top of the Cu atoms when a monolayer of Cu has been 
deposited on the Au(lll) surface [Q. This interpretation is supported by the measurements 
of Shi and Lipkowski [^] and by our own AES results. The total change can therefore not 
be explained by the specified charge-transfer reactions, which would predict a much larger 
charge transfer than the one measured. Shi and Lipkowski [^8|,^ derived an equilibrium 
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electro capillary equation and determined the relative interfacial tension from measurements 
of the electrode charge density, using the back-integration technique. The Cu coverage is 
then simply obtained as the derivative of the relative interfacial tension with respect to 
the electrochemical potential of Cu. They compared the corresponding charge with the 
calculated Cu coverage and concluded that the effective electrovalence of copper varies from 
1.6 to 1.8 in the double-peak range of the electrode potential. These results are consistent 
with those obtained in our present work. 

The total charge transfer measured from CV results varies among the reported experi- 
ments. In the potential range from 120 to 420 mV vs. Ag/AgCl, which includes both CV 
peaks, some of the available coulometric results give values about —350 /iC/cm^ [^^8| -|5T[] 
(also the result of the present work), —460 /iC/cm^ [^], and —470 fiC/cm^ [|^, compared 



with the theoretical value of —444 /xC/cm^ for a monolayer of ions transferring two electrons 
per ion to the Au(lll) surface. According to Refs. ||48|J49[] , variation of the sulfate concen- 
tration from 0.01 mM to 5 mM changes the total charge transfer only by about 5% at the 
negative end of the UPD potential range. Likewise, variation of the Cu(C104)2 concentration 
from 0.01 mM to 5 mM also did not yield significant changes in the total charge. Thus, it 
is highly unlikely that the large discrepancies between the experimentally measured charge 
transfers are due to the differences between the electrolyte compositions used. Rather, we 



believe that the large charge transfers observed in Refs. [^,0] were due to unexpectedly 
large roughness of the electrodes. 

The only systematic studies of kinetic aspects of the UPD in this system of which we 
are aware, are by Holzle, et al. p3| , |69| , |70| . They considered the dependence on the potential 



scan rate of the the separation between the positions of the positive- and negative-going CV 
current peaks (z.e., the hysteresis) and the corresponding maximum current densities. In 
both cases they found a linear dependence on the square root of the scan rate. Additional 
potential-step measurements, in which the potential was changed discontinuously from one 
side of a peak to the other, gave indications of first-order phase transitions and phase change 
via a heterogeneous nucleation mechanism. 



III. EXPERIMENTAL PROCEDURES AND RESULTS 



In this section we detail the experimental procedures employed in the present study, and 
we present our experimental results. In later sections these will be used to formulate our 
lattice-gas model and determine its parameters. 

The combined ultra-high vacuum (UHV) and electrochemistry instrument has been de- 
scribed elsewhere []B7| , [7T| , |7^ . Prior to the electrochemical experiments, the (111) plane of 
the gold single crystal (Johnson Matthey, 99.995% purity) was cleaned by 340 eV Ar^ ion 
sputtering and annealed at around 650 K. The meniscus position of the working electrode 
ensured that only the chosen plane of the gold crystal was exposed to the electrolyte. Fol- 
lowing removal of the electrode from solution, a drop of electrolyte remained hanging from 



the surface. The drop was removed using Teflon tubing terminated with a syringe The 
electrode was next transferred to UHV for characterization by electron spectroscopies. There 
was no carbon contamination on the emersed surfaces, and gold surfaces showed complete 



water desorption [73 



The electrochemical measurements were conducted using conventional three-electrode 
circuitry and an EG&G PAR 362 potentiostat. The working solutions were made of Milli- 
pore water (18 Mfi-cm), ultra-pure grade H2SO4 (Ultrex from VWR), and ultra-pure grade 
HCIO4 (VWR). Solutions were deaerated and blanketed with nitrogen (Linde, Oxygen Free, 
99.99%). The CUSO4 (99.999%, Aldrich) and Na2S04 (ACS grade. Fisher Scientiflc) chemi- 
cals were used as received. Voltammetric current proflles obtained at a scan rate of 2 mV/s 
are shown in Fig. |l]. The separations between the positive- and negative-going CV peak 
potentials were 24 mV for Peak #1 and 44 mV for Peak #2, indicating that kinetic effects 
||53| , |69|j70| may contribute to the shape of the voltammogram to a nonnegligible degree, even 
at this relatively low scan rate. 

The AES was carried out in differentiated mode with 2 eV modulation amplitude, using 
a Perkin Elmer PIII-10-155 cylindrical mirror analyzer (CMA). The spectra were obtained 
using a digital data acquisition system, and smoothed one time using a simple 11-point 



8 



averaging technique. The analysis was performed within 5-8 seconds at a given point on the 
sample, after which the location on the sample was changed for a subsequent 5-8 seconds' 
measurement (to reduce the electron-beam damage). These AES spectra are shown in Fig. |^. 

Quantitative AES analyses of sulfate were performed using a standardization technique 
developed by Wieckowski and co-workers p7||n| , [72| ] . In this procedure, a thick deposit of 
sodium sulfate was formed after water evaporation from gold emersed from aqueous 2 M 
Na2S04. The Au(lll) surface covered by Na2S04 was transferred into UHV and used as 
a standard for the AES characterization of the electrode covered by sulfate up to 1 ML 
coverage. The thick, bulk-like layer of Cu deposited from a solution containing 2 M Cu^"*" 
and 1 mM H2SO4 at a negative potential served as the reference standard for the Cu intensity. 
We used a sensitivity factor approach within a homogenous model of the adsorbed layer, 
resulting from nonlinear equations discussed in Ref. 1^^. The resulting coverages of sulfate 
[obtained from the S(LMM) Auger electron transitions] and Cu [from the Cu(LMM) and 
Cu(MNN) transitions] are shown vs. the electrode potential, together with our simulated 
coverage results, in Fig. |^. 

Our measurements were conducted in a very diluted sulfuric acid solution (0.1 mM), 
containing 1 mM of copper sulfate. The low bulk concentration of copper cations ensures 
that practically none of the Cu-species remains in the emersed film, that could be added 
to the adsorbate coverage measured by AES in UHV. Therefore, no rinsing of the electrode 
surface before emersion to vacuum was necessary. The spot width of the hexagonal (111) 
LEED pattern of the surface indicates a (1x23) reconstruction of the Au(lll) surface, which 
is confirmed by the sharp CV spike in pure H2SO4 solution, seen in Fig. ||(a) at approximately 
440 mV vs. Ag/AgCl [^^. A representative LEED pattern is shown in Fig. 0. 

All measurements were conducted at room temperature. Unless otherwise indicated, the 
electrode potentials are referred to Ag/AgCl. 
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IV. MODEL AND THEORETICAL METHODS 



Our lattice-gas model for the UPD of Cu on Au(lll) in sulfate-containing electrolyte is 
a refinement of the HB model |jl2|,|13,21-Q. It is based on the assumption that the sulfate 



coordinates the triangular Au(lll) surface (STM indicates that the (1x23) reconstruction 
observed by LEED (see Sec. pi]) does not occur in situ Q]) through three of its oxygen 
atoms, with the fourth S-0 bond pointing away from the surface [P!^PT| , P3| , pi ,|5U , oU| (as is 
also the most likely adsorption geometry on Rh(lll) [|T^). This adsorption geometry gives 
the sulfate a "footprint" in the shape of an equilateral triangle with a 0-0 distance of 
2.4 A reasonably matching the lattice constant for the triangular Au(lll) unit cell, 



2.88 A |7g|. The Cu is assumed to compete for the same adsorption sites as the sulfate. 
The configuration energies of the coadsorbed particles are given by the standard three-state 
lattice-gas Hamiltonian (energy function) used, e.g., in Refs. [p3l -|37||: 

(0 (0 (0 



T-Clg = E 

I 

A i i 

Here nf g{0,1} is the local occupation variable for species X [X=S (sulfate) or C (copper)], 
and the third adsorption state ("empty" or "solvated") corresponds to nf=nf=0. The sums 
J^fij) and J2i run over all /th-neighbor bonds and over all adsorption sites, respectively, $xy 
denotes the effective XY pair interaction corresponding to an /th-neighbor bond, and J2i runs 
over the interaction ranges. The term — $gs J^A^^i^j'^^'k expresses three-particle interactions 
between sulfates, involving all second-neighbor equilateral triangles |77|- fr^] . Those bonds 
that correspond to potentially nonzero lateral interactions in the present implementation 
of our model are shown in Fig. ^. For large separations the interactions are assumed to 
vanish, and ^gs* is an infinite repulsion corresponding to nearest-neighbor sulfate-sulfate 
exclusion ("hard hexagons" pT| , p3| , p4| , pO| ) . The change in electrochemical potential when 
one X particle is removed from the bulk solution and adsorbed on the surface is —fix- The 
sign convention is such that $xy>0 denotes an effective attraction, and /ix>0 denotes a 
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tendency for adsorption in the absence of lateral interactions. We emphasize that the $xy 
are effective interactions through several channels, including electron-, phonon-, electrostatic, 
and fluid- mediated mechanisms (see Refs. [p0| , |36| for discussion and references). 

The electrochemical potentials in Eq. (|l]) are (in the weak-solution approximation and 
here given in molar units) related to the bulk concentrations [X] and the electrode potential 
E as 

fc, = ^^l + RT\n^-zsFE 



[S]c 

fic = fi'c + RT\n^-zcFE. (2) 



Here and zc are the effective electrovalences of sulfate and copper, respectively, R is 
the molar gas constant, T is the absolute temperature, and F is the Faraday constant. 
The quantities superscripted with a are reference values which contain the local binding 
energies to the surface. They are generally temperature dependent due, among other effects, 
to rotational and vibrational modes. 

The coverages of sulfate and Cu, 9s and 9c, are the thermodynamic densities conjugate 
to /is and flc, respectively. They are defined as 

9x = iV-iE^f> (3) 

where N is the total number of unit cells in the surface lattice. However, it is experimentally 
observed that the sulfate species remains partially adsorbed on top of the complete monolayer 
of Cu in the negative-potential region, rather than becoming desorbed into the solution 
|^,5B|, as was mentioned in Sec. |T[ Rather than introducing a full, two-layer lattice-gas 



model, we approximate the sulfate coverage in this second layer as 

9^') = a9c(l/3 - 9s) , (4) 

which allows the difference between the first-layer coverage 9s and its saturation value of 
1/3 to be transferred to the top of the Cu layer. The factor a is a phenomenological fitting 
parameter expected to be between zero and one. The total sulfate coverage measured, e.g., in 
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chronocoulometric |^9| and AES experiments, is Gs + 63 , and the total charge transported 
across the surface per unit cell during the adsorption/desorption process thus becomes 

q = -e [zcQc + zs {Qs + QS^)' 
= -e [{zc + azs/3)Qc - zsaQc^s + zsQs)] , (5) 
where e is the elementary charge unit. In the absence of diffusion and double- layer effects 



and in the limit that the potential scan rate dE/dt^O the voltametric current i per 
unit cell of the surface is the time derivative of q. Using differentiation by parts involving 
the relations between the electrode potential and the electrochemical potentials, Eq. (Q), as 
well as the Maxwell relation dQs/dfic = dQc/dfisy we find the current density i in terms of 
the lattice-gas response functions dQx/dfiy: 

dQc 



eF(zl{l-aQc)^ 



+ zcizc - 2azs0s/3) 



+ zsi2zc + azs{l/3 - 63) - azcQc] 



dE_ 
'dt 



(6) 



which reduces to its standard form for a = ]TP| , |5P| -P[] . We emphasize that the derivation 
of Eq. (^) is based on equilibrium statistical mechanics. Very close to a phase transition 
it may require exceedingly slow scan rates to be valid, due to critical slowing-down in the 



case of second-order transitions ||30|-p2| or hysteresis in the case of first-order transitions 
5|,|6|J7|,|8| . 



The following steps are involved in adjusting the lattice- gas parameters to fit the exper- 
imental data. 

A. Ground-state calculations [p^^-p9| to determine parameter ranges that yield ground states 
which correspond to the experimentally observed phases at room temperature . 

B. MC simulations pO|-p^ at room temperature to further adjust the effective interactions 



to obtain reasonable agreement with the shapes of the observed adsorption isotherms and 
CV profiles. 

C. With the effective interactions obtained in step B, determination of the effective elec- 
trovalences so as to give agreement between the dependences of the CV peak positions on 
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the electrol5d;e composition in experiments and simulations. 

D. Iteration of steps B and C as necessary to improve the agreement between the experi- 
ments and simulations. 



A. Ground-State Calculations 

Ordered phases with unit cells up to (3 x 3) were identified by applying the group- 
theoretical arguments of Landau and Lifshitz [^7|-^ . The ground-state energy of each phase 
depends on /is and /2c and the lateral interactions. Among phases with identical coverages of 
Cu and sulfate, those were excluded that have a higher ground-state energy than others with 
the same coverages for all chemically reasonable choices of interaction energies. Out of the 
44 remaining phases, ten were found to be realized as ground states for interactions in the 
range of experimental interest. These are listed in Table where they are identified by the 
notation (XxF)®^. For a given set of interactions, the zero-temperature phase boundaries 
were exactly determined by pairwise equating the energies given in Table |1[ In order to 
easily explore the effects of changing one or more of the interactions, a computer program 
was developed which numerically determines the ground-state diagram by scanning /is and 
fie and determining the phase of minimum ground-state energy p3| . 

The ground-state diagram obtained with the final set of interaction constants used in 
this work is shown in Fig. ^. For large negative fis only Cu adsorption is possible, and 
the phase diagram is that of the lattice-gas model corresponding to the triangular-lattice 



antiferromagnet with next-nearest neighbor ferromagnetic interactions |p4|-[86|]. Similarly, 
in the limit of large positive fis and large negative fie the zero-temperature phase is the 
(v^ X V^)^^ sulfate phase characteristic of the hard-hexagon model [lT2| , pI| , p3| , p4| , |80| . The 



phase diagram for intermediate electrochemical potentials is quite complicated. For /is < 
—22 kJ/mol, no sulfate adsorption occurs in the first adlayer, while for /Zc < — 18 kJ/mol, no 
Cu is adsorbed. The (v^x V^)l^/l mixed-phase region in the upper right-hand part of the 
diagram is relatively large, due to the nearest-neighbor effective attraction between Cu and 
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sulfate, $sc- This is the mixed phase which occurs between the CV peaks. The (v^xv^)o 
phase corresponds to experimental observations in copper-free systems [^,0. This phase is 
enhanced by the fourth-neighbor sulfate-sulfate attraction, $gg , and the {VSy<V^)l^^ ("hard- 
hexagon") phase (which has not been directly observed in this system) is disfavored by ^^Si 
the second-neighbor repulsive trio interactions. The (-\/3xa/7)q^^ phase is also very sensitive 

(2) 

to the next-nearest neighbor interaction $gg . Attractive values of this interaction favor the 
{y/3x\^)l^^ phase, whereas repulsive values favor the (V^xV^)^'^ phase. The experimental 
observations of the (-\/3 x a/T)^^ phase therefore indicate that $gg* must be much weaker 
than the other effective interactions, and we set it equal to zero in our simulations. 



B. Simulations at Room Temperature 

At nonzero temperatures, the thermodynamic and structural properties of a lattice- 
gas model can be obtained by a number of analytical and numerical methods. These in- 



clude mean-field approximations |]87|,|88|, Pade-approximant methods based on liquid the- 



ory |12, 13, 17,21-24], numerical transfer-matrix calculations |33-38|, and MC simulations 
]T|,|2D|,||,|3|-g. The reason for our choice of the numerical MC method is that 

such non-perturbative methods as MC and transfer-matrix calculations provide much more 
accurate results for two-dimensional systems than even quite sophisticated mean-field ap- 
proximations |8^. Yet they are quite easy to program, and with modern computer tech- 



nology their implementation is well within even modest computational resources. (All our 
numerical data were produced on a cluster of IBM RS/6000 workstations.) For models 
with such relatively long-ranged interactions as the present one, MC is better suited than 
numerical transfer-matrix methods. 

The ground-state diagram serves as a guide to the path in the (/xg, fie) plane that the 
room-temperature isotherms should follow. At constant temperature and pB., two factors in- 
fluence the path: the adsorbate concentrations in the bulk, and the effective electrovalences. 
As seen from Eq. (0), /is and /2c depend linearly on E, with proportionality constants — 
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and —zq, which must be determined from experiments. In order to obtain the parts of 
the room-temperature phase diagram of the model that are relevant to the experiments, 
we performed initial isothermal potential scans, using the approximate values zs = —2 and 
zc = +2, which correspond to the full ionic charges of S04~ and Cu^+, respectively. The 
experimental evidence, in particular Refs. HSH^J^, indicates that these 'ideal' values are 



somewhat too large. After obtaining the room-temperature phase boundaries as discussed 
in this Subsection, we therefore re-adjusted our estimates of 2:3 and zq as discussed in Sub- 
ivg below. 



sec. 



The MC simulations used to obtain adsorption isotherms and CV currents at room tem- 
perature {RT = 2.50kJ/mol, corresponding to 28°C) were performed on triangular lattices 
of size Lx L, using a heat-bath algorithm pO|^2| with updates at randomly chosen sites. 
The majority of our simulations were performed with L=30, supplemented with simulations 
for L=45 and 15. Simulation results reported without specifying the system size are for 
L=30. 

In order to avoid getting stuck in metastable configurations (a problem which is exac- 
erbated by the nearest-neighbor sulfate-sulfate exclusion), we simultaneously updated clus- 
ters consisting of two nearest-neighbor sites. Each data point was obtained from a run of 
10^ Monte Carlo steps per site (MCSS), starting from the appropriate ground-state con- 
figuration. Sampling was performed at intervals of 50 MCSS, and the first 4000 MCSS 
were discarded to ensure equilibration. The coverages were obtained as averages over 
these samples, and the nonordering lattice-gas response functions, dQs/dfis, dQc/dfic, and 
dQs/dfic = dQc/dfis, and the adsorbate heat capacity were calculated from the equilibrium 
coverage and energy fluctuations in the standard way [^0| -|32[ . 



The lattice-gas interaction energies were varied until the widths of the CV peaks corre- 
sponding to the room-temperature phase boundaries were in acceptable agreement with the 
experimentally observed peak widths. The corresponding peak positions are superimposed 
on the ground-state diagram in Fig. ^. 
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C. Effective Electrovalences 



Once the room-temperature phase boundaries for the chosen set of interaction constants 
had been calculated, the effective electrovalences could be re-estimated by comparing the 
simulated phase boundaries (CV peak positions) in the (/is,/ic) plane with the dependence 
on the electrolyte composition of experimentally observed CV peak potentials. For this pur- 



concentration of Cu^"*" ions at almost constant bulk sulfate concentration, as summarized in 



A theoretical phase boundary can be represented by a functional relationship, /ig^*^ = 
/ (fiQj. Combining this with Eq. (Q) and differentiating with respect to ln[C], we obtain 



where /' is the derivative of / with respect to its argument, and the subscripts identify the 
particular peak in question. In the parameter region of interest, both the theoretical phase 
boundary and the experimental dependence of the peak potentials on ln[C] are reasonably 
approximated by straight lines. Therefore, combining the results for both CV peaks, one 
obtains two linear equations for the two unknown quantities, zq and zs- Using our simulated 
data, (/{) ^ = — 0.52±0.03 and (/j) ^ = +0.66±0.07, together with the experimentally ob- 
tained (dEPVdln[C])"^ = +42±3 V^^ and (dEf^VdlnlC])"^ = +94±6 y~\ we obtain 
zc = -|-1.68±0.09 and zs = — 1.14±0.16, where the uncertainties only include statistical 
errors. (If, instead, we had used the theoretical ground-state values, {f[)^^ = —0.5 and 
(/2)~^ — +1-0, obtained by equating the ground-state energies in Table. fT|, we would have 
obtained zc = +1.54±0.07 and zs = -0.90±0.12.) 

These values of zc and 2:3 should be considered as averages of weakly potential- and 
concentration-dependent electrovalences in the region of experimental interest. To within the 
error bars they agree with the values recently obtained by Shi et al. |^J51[| by electrocapillary 
techniques. 



pose we used data from Omar et al. who measured the CV peak positions vs. the bulk 



Table |TT|. 
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D. Further Parameter Refinement 

In principle, Steps B and C above can be iterated to further improve the agreement 
between theory and experiment and correspondingly obtain improved parameter estimates. 
However, as shown in the following sections our simulated and experimentally measured 
data agree to a quite satisfactory degree, using the lattice-gas parameters obtained after one 
iteration. Since the iteration process is expensive, both in terms of human and computer 
time, we therefore terminated the loop after the first pass through Step C. Our final MC 
simulations were performed with the effective interactions listed in Fig. ^, using zc=+^-70 
and 2:3=— 1.15. 

V. COMPARISON OF SIMULATED AND EXPERIMENTAL RESULTS 

The room-temperature potential-scan path corresponding to the electrolyte composition 
used in our experiments, was chosen such that the theoretical CV peak separation along 
the scan equals the peak separation in the experimental data shown in Figs. |l](b) and ^(a). 
This scan is indicated by the dotted line labeled "1" in Fig. ^ The simulated CV current 
density [Eq. (^)] along this potential scan is shown together with our experimentally mea- 
sured current (positive-going scan) in Fig. §|(a), and the corresponding integrated charge 
densities [Eq. (^)] are shown in Fig. §|(b). As we pointed out in Sec. Q, there is previous 
experimental evidence that at the negative end of the UPD potential range, sulfate adsorbs 
as a submonolayer on top of the monolayer of Cu, with a coverage e^^^ ^0.2 H. These 
results agree fully with our AES results, which are shown in Fig. ^ The corresponding value 
of a = 0.6 in Eq. (^) was used to obtain the CV current and surface- charge densities shown 
in Fig. ^ and the sulfate coverage shown in Fig. |^ from the lattice-gas simulations. 

With the aid of the ground-state diagram in Fig. || it is easy to analyse the simulated 
and experimental results shown in Figs. ^ and 0. Starting from the negative end, we scan 
in the direction of positive electrode potential {i.e., from upper left to lower right in Fig. 
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Near the CV peak at approximately 200 mV (Peak #2), sulfate begins to compete with 
Cu for the gold surface sites, resulting in a third of the Cu desorbing and being replaced 
by sulfate over a potential range of about 30 mV. Due to the strong effective attraction 
between the sulfate and Cu adparticles, ^sc^ ^ mixed {V^x phase is formed, which 

extends through the entire potential region between the two CV peaks. The (v^x V7)\^/l 
phase which occurs in a narrow band in the ground-state diagram between this phase and 
the (1x1)5* phase that corresponds to the full Cu monolayer, is not sufficiently stable to be 
observed at room temperature. 

As the CV peak at approximately 300 mV (Peak #1) is reached, most of the copper 
is desorbed within a potential range of 10 mV. As it is thus deprived of the stabilizing 
influence of the coadsorbed Cu, the sulfate is partly desorbed, reducing 6s from 1/3 to 
approximately 0.05. The (-\/3x-\/7)o^^ phase in the potential region near 420 mV is consistent 
with experimental observations on copper- free systems [^^. Eventually, more positive 
electrode potentials cause the sulfate to form its saturated (v^xv^)^^ hard-hexagon phase 
near 500 mV. Since this phase has not been experimentally observed in in situ experiments 
for this system, this is probably an artifact of the model at very positive electrode potentials. 

The scenario described here corresponds closely to that proposed by HB 



and the structure proposed for the mixed {^/3x\^)yl^ phase between the two CV peaks 
has recently been strongly supported by Toney et al. with in situ surface X-ray scattering 
||60|| . The nearly linear dependence of the two coverages on each other for Be between 
approximately 0.1 and 0.6, together with the broad maximum of the differential coadsorption 
ratio, [d(6s + 0s )/d6c is close to 0.5 in this region] are shown in Fig. |[ (Also compare 
Fig. P(a) with the insert in Fig. 4(b) of Ref. [^.) These results show that this system 
provides an example of the enhanced-adsorption phenomenon described by Rikvold and 
Deakin |^ . 



Very close to the phase transitions, the CV current profile may be significantly distorted 
from the equilibrium result represented by Eq. (|^) due to kinetic effects. In addition, the 
transition may be rounded by defects or impurities on the adsorbent surface, which can 



have significant effects even at very low concentrations 0]. As a result of these effects, 



the simulated quantities that are most reliably compared with experimental results are the 
CV currents some distance away from their maxima, and the corresponding charge den- 
sities. With this caveat, we find the agreement between our simulations and experiments 
quite good. The shape of the foot region of Peak ^^2 is well reproduced [Fig. ^(a)], and 
the simulated charge transfer across that peak agrees with the experimental value to within 
about 15% [Fig. ||(b)]. Likewise, good agreement is found with the Cu and sulfate coverages 
measured by AES [Fig. |^ . The background current in the ordered-phase region between the 
peaks also shows excellent agreement [Fig. ^(a)]. Somewhat less satisfactory agreement is 
obtained for Peak #1. The simulated peak is somewhat too narrow at the base, and the foot 
on the positive-potential side of the peak, although visible, is lower than its experimental 
counterpart. We find these discrepancies between the simulated and observed current densi- 
ties (which were recorded at a potential scan rate of 2 mV/s and show noticeable hysteresis) 
to be within what one should expect from our essentially equilibrium calculation, even if 
the effective lattice-gas interaction were highly accurate. We therefore believe that a study 
of the effects of kinetics and defects in our model would be very valuable. 

To investigate the effects of changing the electrolyte concentration in the lattice-gas 
model, two additional isothermal potential-scan paths were simulated. We assumed the 
same bulk sulfate concentration as for Path 1 and changed /2c [Eq. @] to simulate bulk 
copper concentrations of 5 mM (Path 2) and 0.2 mM (Path 3), respectively. All three 
paths are shown in Fig. |^ as dotted straight lines. Simulations of current densities and 
coverages corresponding to Paths 2 and 3 are shown in Figs, ^(a) and |^(b), respectively. 
Along Path 2, which corresponds to the highest sulfate concentration, a wider potential 
interval is required to traverse the mixed-phase region; hence the CV peak separation is 
larger than for Path 1. Path 3, conversely, corresponds to the lowest sulfate concentration. 
It exits quickly from the {y/3xy/3)y'^ phase into the disordered low-coverage phase, giving a 
CV peak separation that is smaller than for the other two paths, as well as lower coverages 
of both Cu and sulfate immediately to the positive-potential side of Peak #1. We also 
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note that the experimentally observed asymmetric shape of the foot of Peak #1 is enhanced 
for lower bulk copper concentrations, as is particularly evident from the CV current along 
Path 3, shown in Fig. |^(a). 

VI. FINITE-SIZE SCALING ANALYSIS 

The model studied here is equivalent to a generalized triangular-lattice Blume-Emery- 
Griffiths (BEG) model with interactions of longer than nearest-neighbor range. The 
lattice-gas Hamiltonian in Eq. (|l]) can be transformed into that of the equivalent BEG 
model by the following transformation, 

M = Oc - Os (8a) 

g = ec + Os (8b) 

H =^ific- fis) (8c) 
D = -^{fxc + f^s) , (8d) 

where M is the density conjugate to H, and Q is the density conjugate to —D. 

By comparing the topology of the ground-state diagram shown in Fig. |^ with those for 
the generalized BEG models illustrated in Refs. 0-^,^,0, one can deduce the order and 
universality class of the phase transitions that correspond to CV Peak #1 and #2. 

The transition at Peak #1 corresponds to the phase transition seen in the BEG models 
as D is increased for small nonzero values of H and T. This is a first-order transition in 
which the total coverage Q drops discontinuously from near one to near zero. As is generally 
the case for first-order transitions [^, the adsorbate heat capacity C and the nonordering 
response functions, dM/dH and dQ/dD, diverge proportionally with the system area, L^. 
This divergence is an expression of the discontinuities in M and Q (or equivalently in 0c 
and Gg) and the latent heat associated with a first-order transition. 

The transition at Peak ^2 corresponds to the phase transition seen in the BEG models 
as H is increased for low T and negative D. This is a second-order transition which for a 
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system on a triangular lattice is in the same universality class as the hard-hexagon model, 
namely that of the three-state Potts model [^]. The standard critical exponents in this 
universality class are = § for the correlation length, a = | for the heat capacity, and 
/3 = I for the order parameter [^]. At a second-order phase transition the heat capacity 
diverges with the system size as L"-'^ [pl|], which for this universality class becomes 
The nonordering response functions are coupled to the heat capacity and are expected to 
share the same divergence . 

We have checked these predictions by MC simulations for systems with L=15, 30, and 45. 
The results are shown for Peak #1 vs. Li^ in Fig. |lO|(a) and for Peak #2 vs. L^/s in Fig. |ig(b). 
(To reduce the effects of the finite-size rounding of the phase transitions, we show only the 
scaling results for the response functions corresponding to a direction pointing away from 
the transition line in the phase diagram, —RT{dQ/dD) for Peak ^1 and RT{dM/dH) for 
Peak 7^2. In both cases the response function not shown has the same scaling behavior 
as the one shown, but its magnitude is much smaller.) The results shown in Fig. |10| are 
consistent with the scaling relations predicted above. 



VII. DISCUSSION AND CONCLUSIONS 

In this study we have combined a computational and theoretical lattice-gas modeling 
approach with electrochemical and UHV experiments to study the UPD of Cu on Au(lll) 
single-crystal electrodes in the presence of sulfate. In agreement with other experimental 
and theoretical studies we find that in the potential region between the two sharp CV 
peaks, the electrode is covered by a mixed adlayer of (-\/3x-\/3)i?30° symmetry, consisting 
of 2/3 ML Cu and 1/3 ML sulfate. This ordered-phase region is limited on the positive- 
potential side by a first-order phase transition to a disordered low-coverage phase, followed 
at still higher potentials by transitions to pure sulfate phases. On the negative-potential side 
the mixed phase terminates at a second-order phase transition to a full monolayer of Cu. 
The mixed phase is stabilized by relatively strong, effective attractive interactions between 
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nearest-neighbor copper and sulfate adparticles. 

The lattice-gas parameters were obtained by a three-step procedure. First, a group- 
theoretical ground-state calculation was used to determine the intervals of interaction ener- 
gies which produce ordered phases that agree with the experimental observations. Second, 
MC simulations were used to produce adsorbate coverages and CV current profiles at room 
temperature, and the interaction energies were adjusted until optimal agreement with all the 
available experimental data was obtained. Third, the effective electrovalences were adjusted 
to obtain agreement between the simulations and experiments regarding the dependences 
of the CV peak potentials on the bulk copper concentration. Our final estimates for the 
effective lateral interactions are given in Fig. ^. It is difficult to attach error estimates to 
these values, and there is no guarantee that even better overall agreement could not be 
obtained with a different parameter set. Seen as a formal optimization problem, the param- 



eter adjustment we have performed here by the procedure detailed in Subsecs. [IV A| -D, is 
a daunting task. We are currently studying strategies that could, at least to some extent, 
improve the reliability of the parameter optimization procedure with a manageable computa- 
tional effort However, no first-principles method exists that could calculate the effective 
lattice-gas parameters in this complicated, metal-water-multi-ion system with the degree 
of accuracy needed to calculate thermodynamic and structural quantities for quantitative 
comparison with experiments p7|p8| . We are of the opinion that our approach represents 
the only method for determining effective interaction energies in complicated electrochemical 
systems, which is practically feasible at the present stage of theoretical and computational 
development. 

Despite the cautionary note sounded in the preceding paragraph, we point out that 
our estimated interactions appear entirely sensible from a chemical point of view. Two 
aspects of the interactions are particularly interesting. First, the effective nearest-neighbor 
sulfate-Cu attraction is relatively strong: <l>gc =+4.0 kJ/mol. This is consistent with the 
short Cu-0 bond length (2.15 A) found in recent in situ X-ray diffraction measurements 
]. Second, and perhaps more surprisingly: the nearest-neighbor Cu-Cu interaction is 
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moderately repulsive: $cc =~l-8 kJ/mol. However, interaction energies on surfaces are 
usually much weaker than the interactions between the same atoms in the bulk solid, and 
the presence of the fluid and residual charges on the adsorbed particles may further modify 
the effective interactions to turn attraction into repulsion [|l|. For the present system, the 
conclusion that nearest-neighbor Cu adparticles interact repulsively with each other, was 
already reached by Magnussen et al. 0. We consider it a measure of the robustness of 
our parameter-estimation procedure that we arrived at the repulsive nature of this effective 
interaction independently and only after lengthy attempts at obtaining a reasonable shape 
for Peak ^2 with an attractive nearest-neighbor interaction had failed. 

The residual disagreement between our simulated and experimental results near the CV 
current peaks most likely arise from defects on the electrode surface, and from the use of 
equilibrium statistical mechanics to study in detail a phenomenon that has significant kinetic 
aspects. We believe that the effects of defects and kinetics in our model present promising 
research topics, and we hope to return to their study in the future. 
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TABLES 

TABLE L Representative experimental conditions employed in the present paper and in some 
other voltammetric studies. The peak separations are calculated as averages over positive- and 
negative-going potential scans and are taken from figures or tables in the cited references. 



Reference 



Electrolyte 



Peak separation Scan rate 



This work 
Kolb et al. [|2 
Schultze et al. ^ 
Shi et al. M]M 



Omar et al. [46| 



1 mM CUSO4 0.1 mM H2SO4 

1 mM CUSO4 + 50 mM H2SO4 

1 M HCIO4 + 1 mM Cu2+ + 1 mM S0|" 

1 mM Cu(C104)2 + 5 mM K2SO4 

0.1 M HCIO4 
5 mM CUSO4 + 90 mM H2SO4 



100 mV 
160 mV 
130 mV 
140 mV 

150 mV 



2 mV/s 
1 mV/s 
20 mV/s 
5 mV/s 

5 mV/s 
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TABLE IL List of ordered phases {XxY)q^ and their ground-state energies per site. Only those 
phases are included that are realized as ground states for interactions in the region of experimental 
interest. 





Phase 


Ground— state energy per site 


1 


(1x1)0 





2 


(1x1)? 




3 




-ills + 2/ic)/3 - - 2$g) - 24^(5 - $g - 2$g/3 


4 




-(/is + /ic)/3 - - 4c - 4s - 24s/3 


5 


(a/3xa/3)0/3 


-2/ic/3 - - 2$g) 


6 




-/is/3 - $g - 2$g/3 


7 


(V3xV3)?/3 


-/ic/3-4c 


8 




-(/Zs + 4/ic + $g + 2$g + 9$g^ + 10$g^ + 64g)/5 


9 




-(/is + 3$g)/7 


10 




-(/is + 4s+2$g)/5 
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TABLE III. CV peak positions and separations for different electrolyte compositions at a 



potential scan rate of 10 mV/s, from Ref. [46|. The compositions used were x mM CUSO4 + 
90 mM H2SO4. The peak positions given in the table are in mV vs. the Normal Hydrogen 



Electrode. They are averages over the positive- and negative-going scans in Fig. 4 of Ref. |46| 



taken from Table 2 of Ref. |23|. These results are used in calculating the effective electrovalences 



in Subsec. [V C . 



Peak #1 



Peak #2 



Peak separation 



0.5 
5.0 

50 



444 mV 
496 mV 
554 mV 



303 mV 
327 mV 
352 mV 



141 mV 
169 mV 
202 mV 
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FIGURES 

FIG. 1. (a) Our experimental CV current densities for a Au(lll) electrode in 0.1 mM H2SO4 
at a potential scan rate of 10 mV/s (dotted curve) and in 1 mM CUSO4 + 0.1 mM H2SO4 at 2 mV/s 
(solid and dash-dotted curves). The differences between the latter two curves provide a measure of 
the experimental uncertainties. The arrows indicate the potential-scan directions. Note that the 
current density is proportional to the scan rate. Thus, the copper-free voltammogram corresponds 
to a very small charge transfer per unit potential change, compared to the copper-containing case, 
(b) Integrated charge density for the positive-going scan with Cu^'*"-containing electrolyte (solid 
curve). [The corresponding current density is the dash-dotted curve in panel (a), which for clarity 
is also shown against the right-hand vertical axis in this panel.] The charge density is defined as 
zero at -1-366 mV vs. Ag/AgCl. 

FIG. 2. Auger electron spectra at 3 keV of the Au(lll) electrode covered by uderpotentially 
deposited Cu and coadsorbed sulfate. The emersions were carried out at several electrode poten- 
tials, as indicated vs. Ag/AgCl. These data were used to calculate the Cu and sulfate coverages 
shown in Figs. |^ and ^. 

FIG. 3. The (^/3x^/3)R30° LEED pattern at 91.1 eV electron energy, obtained from the 
Au(lll) electrode emersed from the 1 mM CUSO4 -|- 0.1 mM II2SO4 solution at E=220 mV vs. 
Ag/AgCl, which is in the potential region between the two CV peaks shown in Fig. [l|. 

FIG. 4. The relative positions of Cu (•) and sulfate (A) corresponding to the effective 
interactions in Eq. (Q). The number underneath each bond representation is the corresponding 
effective interaction energy obtained in this work, given in kJ/mol. In particular, the configuration 
corresponding to 'I'ss'' ™ which the corners of the sulfate triangles touch, is forbidden; hence 
the infinitely strong "hard-hexagon" repulsion. The interactions are invariant under symmetry 
operations on the lattice. 
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FIG. 5. The ground-state diagram obtained with the effective interaction constants given 
m Fig. |. The zero-temperature phase boundaries are represented by the sohd hues, and the 
phases are indicated as {X xY)'^^. Superimposed on the ground-state diagram are the positions 
of CV Peak ^1 (fihed diamonds) and Peak #2 (filled squares) at room temperature. The three 
points marked • are low maxima associated with Cu adsorption without coadsorbed sulfate. The 
potential scan path for the simulation corresponding to our experiments with the 1 mM CUSO4 -|- 
0.1 mM H2SO4 electrolyte is shown as the dotted line labeled "1," whereas the lines labeled "2" 
and "3" correspond to additional simulations for 5 mM and 0.2 mM Cu^^ with unchanged sulfate 
concentration, respectively. For all three lines, the upper, left end point corresponds to 120 mV 
vs. Ag/AgCl, and the lower, right end point corresponds to 420 mV. 

FIG. 6. (a) Experimental (dot-dashed) and simulated (solid) CV current densities. The 
experimental curve is the same as the positive-going scan shown in Fig. |l](b). The left-hand 
vertical scale shows the current density for a scan rate of 2 mV/s, whereas the right-hand scale 
shows the current density normalized by the scan rate, in units of elementary charges per mV and 
Au(lll) unit cell, (b) Experimental (dot-dashed) and simulated (solid) integrated charge transfers. 
The experimental curve is the same as shown in Fig. ||(b). The dashed and long-dashed curves 
represent the simulated, partial charge density corresponding to copper and sulfate, respectively. 
The experimental charge density is defined as zero at -1-366 mV. The data points marked x in 
both panel (a) and (b) represent additional simulation results for L=45. They show that finite-size 
effects in the simulations are negligible, except very close to the maxima of the CV peaks. 

FIG. 7. Simulated coverages of Cu, ©c (dashed), first-layer sulfate, @s (dotted), and total 
sulfate, 05-|-G^ (long-dashed), shown together with the corresponding results from our AES 
experiments, □ and O, respectively. 
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FIG. 8. (a) The total sulfate coverage, Gs+Sg , shown vs. Qq- The data points shown as 
crossed error bars represent the AES data from Fig. ^. (b) Finite-difference estimate of the differen- 

(2) 

tial coadsorption ratio d(0s+Og )/d0c as a function of 0c. The ratio exhibits a broad maximum, 
corresponding to the linear section in panel (a), indicative of mutually enhanced adsorption. 

FIG. 9. Simulation results for Path 2 (5 mM Cu^+) and Path 3 (0.2 mM Cu^+) in Fig. |. 
(a) CV current densities for Path 2 (dashed) and Path 3 (solid), (b) Coverages for Path 2 (0c: 
dashed; Gs + Og^'': dot-dashed) and Path 3 (Gc: solid; Gg + 0^^^ dotted). 

FIG. 10. Finite-size scaling plots of the maximum values of the dimensionless adsorbate heat 
capacity, C/R, and nonordering response functions, —RT{dQ/dD) and RT{dM/dH), for system 
sizes L=15, 30, and 45. The straight lines are unweighted least-squares fits to the data points. See 
details in Sec. 0. (a) At Peak #1, showing C/R {O and sohd line) and -RT{dQ/dD) (□ and 
dashed line) vs. L^. The response function has been divided by 10 for easy visual comparison with 
the heat capacity, (b) At Peak #2, showing C/R (O and solid line) and RT{dM/dH) (Q and 
dot-dashed line) vs. L^/^. 
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